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Abstract 

Two-dimensional problem of evanescent wave scattering by dielec- 
tric or metallic cylinders near the interface between two dielectric me- 
dia is solved numerically by boundary integral equations method. A 
special Green function was proposed to avoid the infinite integration. 
A pattern with a circular and a prolate elliptic cylinders, respectively, 
is suggested to simulate the sample and the probe in near-field optical 
microscopy. The energy flux in the midplane of the probe-cylinder is 
calculated as a function of its position. 

The diffraction limit in optics, known to originate from the wave nature 
of light, gives a striking example of a physical restriction being a target of in- 
creasing efforts to overcome. More than a century ago it was realized that the 
wavelength limits the smallest spot within the electromagnetic energy can be 
localized, as well as the smallest details one can optically resolve are compa- 
rable to the wavelength. However, further study showed that these obstacles, 
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in fact, concern a traveling (homogeneous) electromagnetic wave. Unlike, the 
inhomogeneous (also referred to as evanescent) waves, which can not propa- 
gate far away from their source, open the way to clear up the limitations due 
to diffraction and go towards the optics of tiny objects. For instance, the 
nanosized highly-polarizable (i.e. metal) particles well manage to concen- 
trate the light energy within few nanometer range [lj. A near-field scanning 
optical microscopy (NSOM) was suggested to obtain optical signal from the 
objects at nanoscale (see [21 E] and references therein) using sharp tips; the 
latter serve much like an optical antenna [I], which receives the energy of the 
local field and then transfers it to a detector. Thus, the nanophotonics is, 
basically, an optics of evanescent waves, and, consequently, the fundamen- 
tal optical processes (such as diffraction, interference, scattering) are to be 
reconsidered. 

In past two decades a substantial progress is achieved in nano-optics [HIE]. 
However, a significant methodological deficiency persists even for the plain, 
basic problems, like scattering of the evanescent wave by a body. The trouble 
is that evanescent wave can not be considered in isolation from its source (for 
instance, the interface where the total internal reflection takes place), there- 
fore the source is certainly affected by the scatterer as being located within 
a few wavelengths. In paper [7] the general analytical approach is suggested 
that makes possible to do very effective calculations of the evanescent wave 
scattering on a 2D particle (a cylinder) near a flat boundary. In the present 
work we make the next step and consider the problem of two optically cou- 
pled objects placed into the inhomogeneous wave. Our main goal is to get a 
physical insight into the near-field scanning optical microscopy, which min- 
imally involves two small bodies — the studied object and the probe. We 
believe our work is a useful starting point for analysis of particular NSOM 
schemes which will allow for correct extraction of the near field and struc- 
tural information from NSOM data. Keeping in mind this application to the 
realistic configurations, we should focus our attention on the first-principles 
approaches, avoiding restricting assumptions and approximations. 

We start from the Helmholtz equation 



where A is the Laplace operator, k is the wavevector. Consider domain D 
with permittivity e- lTL and its boundary V = <9D. Denote e out the permittivity 
of exterior of D. The Green theorem can be written inside and outside 





2 



domain T>, respectively 



M(r) =i(^ M '- 9 ^) ds ' +M ° (r) - (3) 

Here r ^ T, !K' = IK(r'), d/dn' is the derivative along the internal normal, 
g(r,r') is a fundamental solution to the inhomogeneous Helmholtz equation 

(A + k 2 )g = 5(x-x')5(y-y'). (4) 

It should be noted that function g(r, r ') is not fully arbitrary, namely, should 
satisfy the radiation condition. Hence, the implicit integral over an infinitely 
remote contour in Eq. [3J reduces to the field in the absence of scatterer D, 
denoted as "Kq. 

Let us consider TM wave for which % is the magnetic field. It has to 
satisfy boundary conditions for the field and its normal derivative 

= 0, (5) 

r 

where square brackets denote the jump, e corresponds to either e in or £ out . 
Conditions (jSJ) mean that the magnetic field is always continuous at T, 
whereas its normal derivative has a jump depending on e in , £ out . 

In order to find the field with the help of the Green theorem (J2J), ([3]) 
we need to know "K and d'K/dn at boundary Y. These two independent 
functions satisfy the following coupled equations at r G Y: 

which are obtained by approaching Y from either inner or outer domains. 
Here the fundamental solution inside and outside D is denoted as g in and 
g out , respectively. 

After solving coupled integral equations flBJ), ([7]) the field in arbitrary 
point can be calculated using (T5]), (J3J). This approach is the base for the 
boundary element method (BEM) [HE]. Its advantage consists in diminish- 
ing the problem dimension. For instance, in two-dimensional geometry the 
method deals with one-dimensional contour Y, and then occurs to be very 
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Figure 1: The geometry of light scattering. 



fast and accurate. The method could be extended to several domains. For 
this purpose we must consider integration ([6]), (JZj) along an unlinked man- 
ifold T with corresponding fundamental solution g m inside each. External 
function g out should satisfy the Sommerfeld radiation condition at infinity, 
i.e. be the diverging spherical or cylindrical wave. 

We consider an evanescent wave, Fig.[TJ The plane running wave 3~C(r, t) = 
IKi nc exp (—iut + zkx • r) goes from the dielectric media E\ to other medium 
with permittivity e 2 ; Q\ is the incident angle between the wavevector ki and 
the normal to the boundary. While it is greater than Oq = arcsin ■sJs~ 2 Je\- i the 
angle of the total internal reflection, only the evanescent wave with coordi- 
nate dependence exp(— ny + ik 2x x) penetrates into medium 2, where 

« = -V^isin 2 ^ -e 2 , k 2x = k lx = —y/e~[ sin6 1 , (8) 

c c 

cu, c are the frequency and speed of light. 

Two independent polarization states of falling wave are possible. We 
consider TM-wave with magnetic field vector perpendicular to the plane of 
incidence. This case is more interesting in view of the plasmon resonances 
study, since electric field vector lies in xy plane where the cylinder has a 
finite size. The solution for TE-wave can be considered in the same way. 
The magnetic field of the wave obeys the Helmholtz equation ([Tj), then the 
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BEM method is applicable. However, the problem arises with the infinite 
integration path along x-axis that is hard for numerical calculation. 

To avoid this difficulty we look for the specific Green function G(x, y; x', y') 
satisfying inhomogeneous equation (J3J) in media 1 and 2. Function G depends 
on difference x — x' only, due to translational symmetry. The boundary con- 
dition at y = is 



[G(x,y;x',y')} 



After the Fourier transformation 



ldG(x,y;x',y') 
e dy 



= 0. (9) 

J y=0 



G(x, y; 0, y') = -L f°° G q (y, y')e^ dq (10) 

Eq. (jlj) is reduced to an ordinary differential equation having exponential 
solutions. Using conditions OH]) we obtain the function at y' > in g-domain 

q f[l+r(g)]e^»', y < 0, 

q 2/i 2 y e -^\y-y'\ +r ( g ) e -M2|y+i/'| ) y > . 

Here r(q) = (£i/i 2 — £2^1)/ (£1^2 + £2/^1) is the Fresnel reflection coefficient 
of p-wave at normal incidence [TU], n\ 2 = q 2 — kf 2 . Carrying out Fourier 
transformation ffTOj) of function (TTTJ) at y > we have two terms 



00— iO 

G l = - J e -^ y - y ' l+iq{x - x,) ^-, (12) 

— oo+iO 
00— iO 

q 2 = _ J c -^\y+y'\+ig(x-x') r ( ( l) d 1 _ ^ 

— 00+io 

and G = G\ + G2, where the sign of square root is given by the rule 

The first term ( TT2|) can be calculated analytically and reduces to the 
Green function in the homogeneous space 

Gi(r;r') = Ifl^ (fc 2 p_), (14) 

where .Hq 1 ^ denotes the Hankel function of the first kind [TTJ, p\ = (x — 
x') 2 + (y ± y') 2 . The second term G2 gives the effect of the reflected image 
source. The amplitude of source at each q is equal to the reflection coefficient 
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r(q). Thus along with the point source at (x',y') we have to consider the 
mirror-image source r(q) at (x', —y'). The total field at the upper half-plane 
is the sum of fields generated by the source and its image at each q. The 
Green function automatically takes into account the multiple scattering. The 
Green function of this type was studied for homogeneous waves: spherical 
acoustic, see [T2J [TS], or cylindrical electromagnetic waves [T3] . 

The asymptotic behavior of the Green function in far field can be found 
by the steepest descent method [15]. The stationary point is go = k 2 \x— x'\/p, 
where p = p_ for G± and p = p + for G 2 . The result is a sum of cylindrical 
waves G oc pZ 1 ^ 2 e' lk2p - _j_ rQ p- 1 / 2 e tk 2 p+ ^ w h ere 



Exsinip- v /e 2 (e 1 - e 2 cos 2 (f) 
£1 sin ip + a/ 82(81 - e 2 cos 2 if) 

is the reflection coefficient at q — qo, ip is the polar angle of observation 
counted out from x-direction. The reflection coefficient is ro = — 1 at if = 
and (s/si — y/e~2)/(y/e~i + y/s~2) at if = ir/2. It turns to zero at ip = 



arctan \/£i/e 2 , i.e. at the Brewster angle. 

Green function ()12l) . f[T3"j) . exploited as the external fundamental solu- 
tion for coupled integral equations ©, allows us to avoid the infinite 
integration along Then we can solve the equations for two contours. 

The algorithm has been tested in the case of one contour and homogeneous 
wave when E\ = e 2 . There are analytical formula for circular cylinder [16J and 
numerical calculations for a cylinder with elliptic cross-section [T7] . Compar- 
ison demonstrates the relative consistency within 10~ 4 for iV = 360 panels 
approximating contour T. 

Without scattering body the electric field vector £ in medium 2 has only 
^/-component. The scatterer produces a small component £ x . and the evanes- 
cent wave is partially converted into diverging one. The corresponding Point- 
ing vector 

S = -^-Re (£ x 5T) (16) 

acquires a nonzero y-component, S y , so the energy flux outgoing from the 
plane interface appears. The flux of scattered wave is calculated at distance 
~ 2A, i.e. in the wave zone, and normalized by the average flux of the 
incoming wave in the first medium Si nc = dKf nc /8ir^E~i. Indicatrix of the 
scattering into the upper half-space is shown in Fig. |2J Hereafter the insets 
show the configuration of scattering bodies. The scattering is minimal in 
normal or longitudinal direction and maximal at some medium angles. It is 
the quadruple contribution due to the field of image source (TT5|) and decay 
of the evanescent wave amplitude with y. The angle ip of the first maximum 



6 



105 ™ 75 




Figure 2: Energy flux through a distant semicircle (arb. units) from two 
equal circles with distance between centers L = 0.21 fim as a function of 
polar angle 0<y9<7ratA = 1.512 /zm, R = 0.1 //m, 5 = 0.01 /xm, 
#i = 7r/4, e 2 = 1,^3 = £4 = —91.5 + 10.32, and different S\ = 2 (solid line), 
2.25 (dot-dahed), 3 (dashed). 

increases with £i. For ex — 2 the first maximum is 23° and for £]_ = 3 it is 
41°. The forward-backward asymmetry of indicatrix is a clear evidence of 
violation of the dipole approximation owing to finite sizes (2kR ~ 0.8). 

The method can be applied also to different cylinders. We choose the 
first cylinder (e- m = £3) with round cross section whereas that of the second 
cylinder (e in = e 4 ) is a prolate ellipse with minor semiaxis a = 0.04 /im and 
axis ratio b/a = 10. The major axis is chosen along y-direction. Fig. [3] shows 
the field near the tip as a function of coordinate x at y — 0.15 /xm, when 
the ellipse is moving along to x axis. The pattern in the near-field domain 
is much more complicated than in the wave zone. We see the slight low- 
frequency oscillations at the right side and deep high-frequency at the left. 
The physical nature of the oscillations consists in the interference between 
falling evanescent wave and the diverging wave scattered by the circle. They 
are counter-propagated at the left and co-propagated at the right. Their 
spatial frequencies are different k± = k x ± k 2 . However, in considered case 
k x = fcisin#i ~ k 2 , then the interference oscillation at the right side has a 
nearly zero frequency. The slight interference pattern at the right in this case 
is caused by the scattering by the tip and these oscillations vanish without 
the tip. 

Fig. H] shows the energy going along the major axis through the ellipse 
middle cross-section as a function of coordinate x, namely, the component 
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Figure 3: The magnetic field at the bottom point of ellipse as a function of 
coordinate x at distance y = 0.15 from the plane (solid line). Parameters are 
A = 1.512 /mi, R = 0.06 /mi, 5 = 0.01 /mi, e x = e 3 = e 4 = 2.25, e 2 = 1, 0i = 
7r/4. The same without the ellipse (e 4 = 1, dashed). 




Figure 4: Energy going through the central cross section in ^-direction of the 
ellipse as a function of its horizontal position x at distance from the plane 
y = 0.15 (solid line), 0.35 (dashed), 0.50 /mi (dot-dash). Parameters are 
A = 1.512, R = 0.06, S = 0.01 /mi, e 1 = e 3 = e 4 = 2.25, e 2 = 1, B x = tt/4. 
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Sy of vector f[T6"j) averaged over the horizontal cross-section. The vertical 
distance of the tip from the interface is d — 0.15 -f- 0.50 /mi, then the min- 
imal distance between the tip of ellipse and the circle starts from 0.02 /mi. 
We see interference oscillations in the coordinate dependence, the amplitude 
decreasing with the distance between the tip and the plane. As follows from 
Fourier transform fill I) , the higher spatial harmonics with q 2 > k\ decay 
exponentially with distance like exp(—fi 2 y)- Then the small details are not 
visible in the far-field pattern. However at ky <C 1 the exponent is not neg- 
ligible, and then the small-scale details become apparent. The near field can 
be observed if one extract the signal and transfer it to the far zone. In our 
calculation the stretched ellipse plays a role of such a transmitter. 

The energy flux through the far central plane is a simplest (two-dimensional 
in our case) model of NSOM [21 El E]- The curves in Fig. @] correspond to 
instrumental function of the microscope. Although, this statement should 
not be taken literally. The multiple scattering leads to back influence of the 
probe to the object then it is not a usual linear function of response. 

The coupled boundary equations describing the scattering of the evanes- 
cent wave are solved for two cylinders. The asymmetry of indicatrix and 
oscillations in the coordinate dependence are observed. The BEM with pro- 
posed Green's function is rather general and applicable for any contour T 
or several contours. It can be extended also for 3D geometry. The Green 
function could find applications in other calculations, e.g. volume integral 
equations including the Born series and discrete dipole approximation. 
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is supported by the Government program NSh-4339.2010.2, program # 21 
of the Russian Academy of Sciences Presidium, and interdisciplinary grant 
#42 from the Siberian Branch of RAS. 
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